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Travelling-wave solutions of the inviscid Burgers equation having smooth initial wave 
profiles of suitable shapes are known to develop shocks (infinite gradients) in finite times. 
Such singular solutions are characterized by energy spectra that scale with the wave number 
k as k~ 2 . In the presence of viscosity v > 0, no shocks can develop, and smooth solutions 
remain so for all times t > 0, eventually decaying to zero as t — > oo. At peak energy 
dissipation, say t = t*, the spectrum of such a smooth solution extends to a finite dissipation 
wave number k v and falls off more rapidly, presumably exponentially, for k > k v . The 
number N of Fourier modes within the so-called inertial range is proportional to k v . This 
represents the number of modes necessary to resolve the dissipation scale and can be thought 
of as the system's number of degrees of freedom. The peak energy dissipation rate e remains 
positive and becomes independent of v in the inviscid limit. 

In this study, we carry out an analysis which verifies the dynamical features described 
above and derive upper bounds for e and N. It is found that e satisfies e < u 2a ~ l \\u* H^, 1 (■ 
where a < 1 and = u(x, £*) is the velocity field at £ = £*. Given e > in the limit v — > 0, 
this implies that the energy spectrum remains no steeper than k~ 2 in that limit. For the 
critical k~ 2 scaling, the bound for e reduces to e < \/3k {{uqW^ ||mo|| 2 , where k Q marks the 
lower end of the inertial range and «o = u(x,0). This implies iV < \/3L ||wo|| 00 /i', where 
L is the domain size, which is shown to coincide with a rigorous estimate for the number 
of degrees of freedom defined in terms of local Lyapunov exponents. We demonstrate both 
analytically and numerically an instance where the k~ 2 scaling is uniquely realizable. The 
numerics also return e and t*, consistent with analytic values derived from the correspond- 
ing limiting weak solution. 



1 Introduction 



In 1948 Burgers 1 introduced the equation 

u t + uu x = vu xx (1) 

as a model for fluid turbulence. Here, u(x,t) is a one-dimensional velocity field and v > 
plays the role of viscosity in a usual fluid. On the one hand, this model captures the 
two most fundamental features of fluid dynamics by its quadratic advection and viscosity 
terms. On the other hand, Eq. ([TJ) lacks a pressure term, thus governing a hypothetical 
compressible fluid without pressure. The absence of a pressure-like term makes Eq. (pQ) 
integrable by the Cole-Hopf method. 2 ' 3 This renders Eq. ([T]) and its generalization to 
higher dimensions poor models for fluid turbulence. Despite this apparent shortcoming, 
the Burgers equation has been widely studied for a variety of applications. 4-13 

The development of shock waves or discontinuities (infinite gradients) from suitable 
smooth initial velocity profiles is an intrinsic property of the inviscid Burgers equation. 
Given a different iable initial profile u(x, 0) = Uq(x), Eq. ([T]) with v = is implicitly solved 
by the travelling-wave solution 

u(x, t) = ito(6) = u o( x — ut). (2) 

By taking the spatial derivative of Eq. ([2]) and solving the resulting equation for u x one 
obtains 

U x = -^-r, (3) 

where u' (£) denotes the derivative of «o(0- ^ follows that u x diverges {u x — > — oo) provided 
that «o(0 < f° r some £• The earliest time t = T for this to occur is T = —l/u' (xo), 
where u' {xq) is the steepest slope of u$(x) occurring at x = xq. This steepest slope travels 
at the speed uq(xq) and gets ever steeper as t — > T, becoming infinitely steep when t = T 



at x = Xq + u (x )T = x — u (x ) /u' Q (x ). In summary, the space-time coordinate of the 
shock is 

Such a singular solution is characterized by an energy spectrum E{k) that scales with the 
wave number k as E(k) oc k~ 2 , which is the spectrum of a step function. 

Under viscous effects, the would-be shock is suppressed, and the solution remains 
smooth and decays to zero in the limit t — > oo. This statement is true however small 
the viscosity. This means that the maximally achievable (peak) energy dissipation rate, 
hereafter denoted by e m , remains positive in the inviscid limit v — > 0. For fixed v > 0, the 
velocity gradient \u x \ can achieve a finite maximum only. Presumably, the corresponding 
energy spectrum would retain the k~ 2 scaling up to a finite dissipation wave number k u , 
around which the dissipation of energy mainly takes place and beyond which a more rapid 
decay, probably exponential decay, occurs. Given this scaling, e m scales as vk v . It follows 
that the number N of Fourier modes within the wave number range k < k v , the so-called 
inertial range, is 

iV oc k v oc — , (5) 

for dimensionally appropriate proportionality constants. This is the number of modes 
necessary to resolve the dissipation scale and can be considered the system's number of 
degrees of freedom. 

In this study, we carry out an analysis that quantitatively confirms the dynamical 
features described above. It is found that e m satisfies 



^ 2a— 1 II i|2(l-a) 



2 



(-A) a /V " , (6) 



where a < 1, A is the Laplace operator, u* = u(x,t*) is the velocity field at the time of 
peak energy dissipation t = t*, and ||-|| and ||-|| denote L°° and L 2 norms, respectively. 



Given that e m > in the limit v — > 0, this result implies that the energy spectrum 
E(k,t*) becomes no steeper than k~ 2 in that limit. For this critical scaling, e m is found to 
satisfy e m < V3k Hw*^ ||w*|| 2 < V^&o ll^olloo H^oll 2 ; where u = u(x, 0) and fc is the wave 
number that marks the lower end of the energy inertial range. This result further implies 
k u < y/3 HmoIIoo / v - It follows that N < y/3L {{uoW^ / is, where L is the domain size, which is 
shown to coincide with a rigorous estimate for the number of degrees of freedom defined in 
terms of local Lyapunov exponents. Note that one can identify the upper bound for iV with 
the Reynolds number Re as in the case of a real fluid. Thus, the system's number of degrees 
of freedom scales linearly with Re. We demonstrate both mathematically and numerically 
an instance where E(k, £*) oc k~ 2 is uniquely realizable. The numerics also return the values 
of e m and t* which are consistent with those derived from the corresponding limiting weak 
solution. 

2 Energy dissipation and dissipaton wave number 

For simplicity, we consider periodic solutions of Eq. (pQ) having period 2ttL and vanishing 
spatial average. The usual L p norm of u (and of its derivatives), for all p > including 
p = oo, is defined by ||w|| = (|w| p ) 1 / p , where (•) denotes a domain average. The advection 
term of the Burgers equation conserves \\u\\ . Under viscous effects, ||u|| decays for p > 1 
and is governed by. 

|iNi P =-^-i)Nt p (i«r 2 ^)- (7) 

Since we are dealing with L 2 and L°° norms only, we omit the subscript p = 2 in the former 
for convenience. The decay of the energy ||w|| 2 /2 is governed by 

1 d II II 2 || 11=2 (o\ 

-— it = — v \\uJ\ . (8) 
2dt 11 11 11 11 v ; 

This section is mainly interested in optimal estimates for the decay rate v \\u x \\ , particularly 

in the limit of small is, and related issues concerning the energy inertial range. 



The governing equation for the velocity gradient u x is 



u xt + uu xx + u n 



(9) 



By multiplying Eq. (jTJ by u xx (or Eq. (|UJ) by u x ) and integrating the resulting equation 
over the domain we obtain the evolution equation for the mean-square velocity gradient 
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where the inequality is straightforward. The final line of Eq. ffTOj) can be used to derive 

1 1 2 

an upper bound for the energy dissipation rate • For this purpose, consider the 

inequality (see Eq. (7) of Ref. 14) 
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where a < 1 is a parameter, which can be varied for an optimal bound, and A is the 
Laplace operator. The fractional derivative (— A) Q / 2 is a positive operator and is defined 
by (— A) a / 2 u = k a u, where (— A) a / 2 u and u are the Fourier transforms of (— A) a ^ 2 u and u, 
respectively. Upon substituting Eq. (TTTi) into Eq. fTTUT) and noting that d \\u x \\ 2 /dt = at 
the time of peak energy dissipation t = t*, we can deduce that 



e < v 2a ~ l \\u \\ 2{1 ~ a) 

c m _ " II lloo 



-Ayi 2 u* 



(12) 



where is bounded by its initial value, but (— A) a ^ 2 u* can be large, depending on 

both E(k,t*) and a. In section IV, we demonstrate both analytically and numerically that 
in the limit — ^ 0, t* is independent of v and, in general, not related to the singularity 
time T of the corresponding inviscid solution. 



Equation (TT21 confirms the fact that e m < oo (and hence \\u x \\ < oo) for v > as one 
can set a = and obtain e m < 11^*11^ ll^*!! 2 /^ < ll M o|lL ll^oll 2 /^- This bound can be highly 
excessive, and a more optimal estimate is possible by varying the "optimization" parameter 

2 

a within the permissible range a < 1. Observe that the spectrum of (— A) a ^ 2 u /2 is 
k 2a E(k). So, if the energy spectrum E(k, £*) is strictly steeper than k~ 2 , then (— A) a ^ 2 u* 
is bounded for some a > 1/2. If this were the case for all u, including the limit v — > 0, 
then the upper bound for e m in Eq. f|T2l would vanish, thereby contradicting the fact that 
e m > in that limit. This rules out energy spectra steeper than k~ 2 . In section IV, we 
mathematically demonstrate an instance where energy spectra shallower than k~ 2 are also 
ruled out. Thus the scaling k~ 2 is uniquely realizable. This suggests that in general, the 
most plausible scenario is that in the inviscid limit, E(k,t*) approaches the k~ 2 critical 
scaling. 

Now, suppose that E(k) = Ck~ 2 /2, for k G [fco,^], where C > is a constant. Note 
that k Q is not necessarily the lowest wave number 1/L. We then have ||w|| 2 = C k~ 2 dk, 
so C = k \\u\\ 2 . Thus, E(k) = k \\u\\ 2 k~ 2 /2. For this case, a direct estimate of the ratio 
ll^x || / ll^xx || is 

IN!! = v^olM 2 . (13) 

W^xx || 

By applying this equation to and substituting the resulting estimate into Eq. ffTUl) we 
deduce the upper bound 

e m < a/3/c IKIIoo IKII 2 - || w Hoo ||mo|| 2 • (14) 

We find later by an example for the parameter values Hwolloo — 1, H M o|| 2 — 1/2 and k = 1 
that e m = 0.1061, which gives us a sense of the sharpness of the derived upper bound 
V^o | [ tio II cc ll M o|| 2 = \/3/2. The dissipation wave number k u , which marks the end of the 
k~ 2 inertial range, is found to satisfy 



K < ^JNL . (15) 



V 



It follows that the number N of Fourier modes within this inertial range is bounded by 

N < ^JNk = Re , (16) 

where Re is the Reynolds number. Note that this estimate also includes the modes cor- 
responding to k < k . The linear dependence of N on Re is interesting and is rigorously 
verified, without reference to E{k ) t 1f ) 1 in what follows. 

3 Lyapunov exponents and number of degrees of free- 
dom 

This section derives a rigorous estimate for the number of degrees of freedom, which is 
defined as the minimum number of greatest local Lyapunov exponents (of a general tra- 
jectory in phase space) whose sum becomes negative. This number, denoted by D, is the 
dimension of the linear space (spanned by the corresponding Lyapunov vectors), which can 
adequately "accommodate" the solution locally, and is essentially the so-called Lyapunov 
or Kaplan- Yorke dimension. 15 ' 16 Its estimate is found to agree with that for iV obtained 
earlier in the preceding section. This agreement is not coincidental and can be considered 
as analytic evidence for the expected k~ 2 energy spectrum used in the estimation of N. Like 
N, D can be thought of as the number of Fourier modes necessary to resolve the steepest 
velocity gradient during the course of evolution, particularly around t — t*. We follow the 
procedure formulated by Tran and Blackbourn 17 in the calculation of the number of degrees 
of freedom for two-dimensional Navier-Stokes turbulence. For a detailed discussion of the 
significance of D, see Refs. 17 and 18 and references therein. 

Given the solution u(x,t) starting from some smooth initial velocity field uq(x), con- 
sider a disturbance v(x,t) satisfying the same conditions as u(x,t), i.e., periodic boundary 
condition and zero spatial average. The linear evolution of v(x,t) is governed by 

v t + uv x + vu x = vv xx . (17) 



The governing equation for the norm \\v\\ is 



Ml ^ IN = —(v(uv x + vu x )) - v \v x f 
= (uvv x ) - v \\v x f 



< I^IL \\v\\ \\v x \\ - u\\v x \ 



< INIL INI \\v x \\ -v\\v x \ 
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where we have used (v 2 u x ) = —2{uvv x ) by integration by parts and the inequalities are 
straightforward. Dividing both sides of Eq. f|T8|) by ||t> || 2 yields 
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where A is the exponential rate of growth (A > 0) or decay (A < 0) of the disturbance norm 



\v\\. 



The set of n greatest local Lyapunov exponents {Ai, A 2 , • • • , A n } and the corresponding 
orthonormal set of n most unstable disturbances {v 1 , v 2 , • • • , v n } can be derived by succes- 
sively maximizing A with respect to all admissible disturbances v subject to the following 
orthogonality constraint. At each step i in the process, the maximizer v is required to 
satisfy both \\v\\ = 1 and (w 3 ') = 0, for j — 1, 2, • • • , i — 1, where f- 7 is the solution obtained 
at the j-th step. Since each normalized solution (Aj,i/) satisfies Eq. f|T9|) . we have 
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where c is a constant independent of the orthonormal set in question. In Eq. (1201) . we 
have applied the Cauchy-Schwarz inequality Yh=i IMII — ( n Sr=i IMI| 2 ) 1//2 an d use d the 



estimate 

which is a consequence of the Rayleigh-Ritz principle. By this principle, the left-hand side 
of Eq. (|2ip is not smaller than the sum of the first (i.e., smallest) n eigenvalues of —A. 
These eigenvalues are l/L 2 ,2 2 /L 2 ,---,n 2 /L 2 and sum up to n{n + 1) (2n + 1) / (6L 2 ) . Hence, 
Eq. (121]) follows with c tending to for large n. Now the condition Yh=\ Aj < is satisfied 
when n > cL ||wo|loo l v - ^ follows that 

D<c^. (22) 

This estimate agrees with the upper bound (fT6l) for N, which was derived by assuming 
the energy spectrum E(k) oc k~ 2 . This agreement provides us with confidence in the 
plausibility of the k~ 2 scaling. 

The term on the right-hand side of Eq. fl2"2l is the Reynolds number Re defined earlier 
with c = v^- Thus D scales linearly with Re. For a comparison, D scales as Re(l + 
In Re) 1,/3 and Re 9,/4 for two-dimensional and three-dimensional turbulence, respectively. The 
former has recently been derived 17 while the latter is a classical result deduced from the 
Kolmogorov theory. These scalings reflect the intrinsic characteristics that the dynamics 
of the two-dimensional vorticity gradient and three-dimensional vorticity are effectively 
linear and quadratically nonlinear, respectively 18,19 The present finding of exactly linear 
dependence of D on Re is somewhat unexpected as the Burgers velocity gradient dynamics 
are quadratically nonlinear, just as the three-dimensional vorticity dynamics. Nonetheless, 
this is not a total surprise if the dimension of the physical space, which plays a significant 
role in the scaling of D with Re, is taken into account. 18 Note that in all three cases, D 
scales linearly with the domain volume, given all else fixed. This is in accord with the 
notion of extensive chaos. 20-22 The linear scaling of D with Re for the Burgers case is fully 
justified in the numerical simulations reported in the next section, where we observe that 



the ratio D/Ke is best kept fixed (at order unity) for various resolutions. Hence, doubling 
the resolution (i.e., doubling D) allows the viscosity to be halved, given all else fixed. This 
allows the exponential dissipation rate vk 2 at the truncation wave number to grow as Re. 
On the other hand, this same linear scaling of D with Re in two-dimensional turbulence 
means that numerical simulations can be performed using a fixed dissipation rate vk 2 at 
the truncation wave number, for different resolutions. Thus, doubling the resolution (i.e., 
quadrupling D) allows the viscosity to be reduced by the factor 1/4. This fact is well known 
to numerical analysts. The scaling of D as Re 9 ^ 4 in three-dimensional turbulence implies 
that the dissipation rate vk 2 at the truncation wave number should be proportional to 
Re 1/2 . This means that doubling the resolution (i.e., octupling D) can allow the viscosity 
to be reduced by the factor 2~ 4 / 3 . 

4 A case study 

In this section we analytically and numerically consider an example that confirms the results 
derived in the preceding sections. In addition, we prove that no power-law energy spectra 
other than k~ 2 are realizable, thus giving an exact result of the slope of E(h,t*) rather 
than a constraint for this particular case. We also determine by numerical simulations 
the viscosity-independent maximum dissipation rate e m and the corresponding time t = £* 
when this occurs. The numerical values of these dynamical parameters agree with those 
derived from the corresponding limiting weak solution. 

4.1 Analytical consideration 

We consider the periodic domain [— ir, ir], i.e., L — 1, and uq(x) = — sin x. This initial profile 
was used in a computational study 13 of the Burgers equation, using 4096 grid points. In 
the next subsection, we report results from simulations using up to 4 x 10 4 Fourier modes. 
It can be readily seen that Eq. ([1]) admits odd functions as solutions. In other words, if 



f(x,t) is a solution, then f(—x,t) is also a solution provided that f(x,t) = —f(—x,t). 
Hence, for the initial profile under consideration, u(x,t) remains odd for all t > 0. We can 
then express u(x,t) in terms of an odd Fourier series: 

u(x, t) — y^-Ufc(t) sinfca:, (23) 
k 

where k — 1, 2, 3, • • • are the wave numbers. The gradient m x is given by 

u x (x, t) — 7^ kuk(t) cos fcx. (24) 
fe 

The origin is "stationary" and has the steepest negative slope, initially equalling —1, which 
is given in terms of Uk by 

u x (0,t) = J2ku k (t). (25) 

k 

The third derivative u xxx (0,t) is 

u xxx {0,t) = -J2k 3 u k (t). (26) 

k 

By substituting Eqs. (1231) and (J23j) into Eq. (jUJ) one obtains 



(27) 



In the inviscid case, u x (0,t) — > — oo as t — *■ T — 1. This can be seen either by solving 
Eq. (1271) with v = or directly from Eq. (j4]). Figure [1] illustrates the viscous solution (for 
v = 0.02) at a few selected times before, near and after the inviscid singularity time (t — 1). 
The evolution of the Fourier coefficients «&(£) is governed by 

—u k = -u 2 k/2 T- ^ u m u e - uk 2 u k , (28) 

m±i=k 

where the sum is over all pairs of wave numbers m and £, including m = £ = k/2 when 
k is even, satisfying the triad condition m ± i = k. Within each individual wave number 
triad, the energy is conservatively transferred from each of the two lower wave numbers 
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Figure 1: A viscous solution to Burgers equation starting from u(x, 0) = — sinx, for v = 
0.02, and shown at times t = 7r/8, 57r/16 ~ 1, 7r/2, 37r/4 and tt. 

to the third and higher wave number or vice versa. It can be seen that all wave numbers 
are initially excited in such a way that Uk < 0. Plausibly, no particular modes would 
become completely depleted of energy during the subsequent evolution. This means that 
Uk does not change sign and remains negative. This fact is verified below in the numerical 
simulations. As a consequence, the transfer of energy to ever-smaller scales is irreversible, 
and each Fourier mode contributes to the steepness of the slope u x (0,t) as there are no 
cancellations in the sum Y^k kuk- The nonlinearity can be said to operate at "full strength," 
without "depletion." This is consistent with the fact that u x (0,t) quickly diverges if v — 0; 
indeed 1^(0, 1) = — oo. This observation prompts us to take Uk < for all k in what follows. 

Consider the inertial range scaling Uk = — c 7 /c -7 , for < 7 < 3/2 and c 7 > 0, which 
corresponds to the energy spectrum E{k) = c 2 k~ 21 /2. By substituting this scaling for Uk 
into the right-hand side of Eq. (|2"T|) we obtain 

<9t V (2-7) 4-7 

_ 1+7 / (3-2 7 )6 m _ c%kl-^ 
- » W4- 7 ) (2-7) 2 

where e m = uc^kl' 21 / (3 — 27) has been calculated from the above spectrum. The fact 
that both J2k k u k — 00 and < e m < 00 as k v — > 00 requires 7 = 1, which is the only 
possibility allowed by Eq. (|2^|) . Indeed, if 7 > 1 (which has already been ruled out in 





Figure 2: A schematic description of energy loss after wave breaking at t = 1 for the 
travelling-wave solution u = — sin(x — ut) of the inviscid Burgers equation. The energy 
dissipation rate is U 3 (t)/ (3tt), where 2U(t) is the shock width. This rate is zero upon wave 
breaking and grows to its maximum of 1/ (3tt) at t = tt/2. 

general) , then the second term in the brackets of Eq. (1291) could be made arbitrarily small 
for sufficiently large k u and the right-hand side would become positive. This contradicts the 
fact that Yjk k u k —oo. On the other hand, if 7 < 1, then the second term in the brackets 
of Eq. (1291) could be made arbitrarily large for sufficiently large k u and the right-hand side 
would become negative. The gradient at the origin J2k ^ u k would diverge for k u < 00, 
which is not possible. 

We now consider the energy dissipation rate in the inviscid case due to the lack of 
smoothness of solution after wave breaking at t = 1. This consideration allows us to 
determine the energy dissipation rate, among other things, of the viscous case in the inviscid 
limit. For t > 1, the travelling- wave solution becomes multivalued in a neighborhood of 
x = as the respective portions u > and u < of u cross over the vertical axis, invading 
the region x > and x < (see Figure [2]). Consider the weak solution consisting of two 



disconnected travelling- wave branches u + (x,t) and u_(x, t) given by 

, s f - sin(:r - u+£) for -tt < x < , . 

u +( x > t ) = { for 0<x<, (30) 

and 

, ,, i — sin(x — U-t) for < x < n 

u -( a; '*) = ( far -*<*<(). (31) 

These terminate on the vertical axis at u + (0, t) = U(t) and u-(0,t) = —U(t), where the 
(half) shock width U(t) is given implicitly by U — sm(Ut). Evidently, lim t _ >1+ U(t) = 
and U(n/2) = 1, the latter of which is the global maximum. The evolution of the energy 
corresponding to this solution is governed by 

1 c\ 1 / fO I-* 

1 " " l|2 _ 1 ..2 /.. \ j„ , / ..2 
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The energy dissipation rate f/ 3 / (37r) tends to zero as £ — > 1 + and achieves its maximum 
of l/(37r) at t = 7r/2 when U(n/2) = 1. For t > 7r/2, this rate decreases monotonically 
to zero as t — > oo. Since the viscous solution approaches this (unique) weak solution in 
the limit v — > 0, the limiting energy dissipation rate for £ > 1 is £/ 3 /(37r). The maximum 
dissipation rate corresponds to U — 1, i.e., e m = l/(37r), occurring at £ = t* = ir/2. Note 
that £* differs from T and is the time for the extrema (initially at x — ±7r/2) to arrive at 
the stationary shock position x — 0. In the next subsection, we recover both values of e m 
and £* with high precision by numerical simulations. 

An interesting feature of the present problem is that in the inviscid limit the energy 
commences its decay from £ = 1 while the maximum velocity does so from £ = ir/2, upon 
which the energy dissipation reaches its peak. This lag in the dissipation of IHI^ can 
be readily appreciated by the following observation. For the energy, the dissipation rate is 



dominated by |w x (0, t)\, which becomes sufficiently large at t = 1, upon which the transition 
between nondissipative and dissipative phases takes place. For the maximum velocity, by 
taking the limit p — > oo of Eq. ([7]) we obtain 

4 \\u\\ = -v lim (p - 1) \\u\\l~ p (\u\ p - 2 ul). (33) 

The dissipation rate on the right-hand side of Eq. (133]) is dominated by \u x \ in the vicinity 
of the maximum velocity. Evidently, as the maximum velocity approaches the vertical axis, 
\u x \ in its vicinity becomes greater (see figure [1]). The transition between inviscid and 
viscous dynamics of WuW^ at t = ir/2 implies that \u x \ in this vicinity is not sufficiently 
large until t = tt/2. A similar behavior has been observed numerically in two-dimensional 
turbulence, whereby the vorticity supremum remains virtually unchanged until (and even 
after) the dissipation rate of the mean square vorticity has achieved its maximum value. 23 

The weak solution provides a convenient way for calculating the dissipation rate d IM)^ / dt 
for t > tt/2. In the limit u ^ 0, one can identify (Hl^ with U = sin(Ut). By taking the 
time derivative of this expression and solving the resulting equation for dU /dt = dWuW^/dt 
we obtain 
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(34) 



In the present example, — u x (x, 0) peaks at an isolated point, namely at x = 0. The 
weak solution is a step function with U(T) = and the energy dissipation rate tends to 
zero as t — ► T + . Similarly, consider a smooth initial profile u(x, 0), for which — u x (x, 0) 
achieves a positive maximum at a finite number, say iVo, of isolated points. Such a profile 
evolves into a piecewise smooth solution having iVo steps, each with U(T) = 0. For this 
case, the energy dissipation rate also tends to zero as t — > T + . When the said maximum 
occurs over an extended interval, say [xi,^], then U(T) = (x\ — x 2 )u x (xi,0) > 0. The 
energy dissipation rate upon wave breaking jumps from zero to a positive value. 








-2- 



l °S, (- u k)- 



v = 7/76000 



-4- 



-6- 







2 



3 



4 



5 



Io 8w k 



Figure 3: Spectra (— w fc vs k) at t = 7r/2 for the three smallest values of viscosity considered, 
i/ = 1/4000, 1/8000 and 1/16000 (computed at resolutions k mSLK = 10000, 20000 and 40000 
respectively). Note, the spectra differ negligibly except in their high wave number tails, 
and are well fit by a k~ x slope in the inertial range. 

4.2 Numerical results 

We now turn to results of a numerical analysis of the Burgers equation. We have simulated 
the initial value problem described by Eq. (1281) . where Wi(0) = —1 and ttfc(O) = for k > 1, 
for several different resolutions up to fc max = 4 x 10 4 . For this given initial condition and 
c = a/3, Eq. (1221) becomes D < y/3/ v. The viscosity v = 2.5//c max has been chosen in accord 
with this estimate to ensure that A; max lies well within the dissipation range. Our choice 
turns out to yield adequate dissipation, thus providing evidence for the sharpness of Eq. 
( |22l) . We have used a standard 4th order Runge-Kutta method with the viscosity exactly 
incorporated through an integrating factor. The adapted time step St = — 0.01/ Y^,k k> u k 
has been used to account for the highly sensitive nature of the problem when t ~ i*. 

Figure [3] shows the plots of log[— versus log/c for the three highest-resolution 
simulations. These exhibit a clear slope of —1 in the inertial range, thus implying the 
scaling k~ 2 for the energy spectrum. Evidently, the inertial range becomes wider for higher 
Re, and a careful inspection of data also shows a clear trend that the inertial range becomes 
shallower, approaching the critical scaling k~ l as expected. 
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Figure 4: Evolution of the energy dissipation rate e(t) for a series of six simulations differing 
in v by factors of 2 (the extreme values of v are indicated). Also, the inviscid singularity 
time it = 1) is indicated by the vertical dashed line. 

Figure |H shows the evolution of the energy dissipation rate e(t) = v \\u x \\ 2 = vJ2k ^ 2-u l/2 
from t = to t = 7r. The dissipation rate remains small for t < 1 (evidently tending to 
zero in the inviscid limit), only to grow considerably when t = 1, consistent with the 
result (1521 for the limiting weak solution. This rate continues to increase for t > 1 and 
achieves a maximum at t — i* = 1.571, which is very close to the analytic value n/2. This 
value of £* has been observed to be very robust with respect to independent variations of 
the Reynolds number and the time step. The maximum dissipation rate is e m = 0.10605 
for the three highest Reynolds numbers. This suggests that the convergence of e m as 
v — > is rapid. Indeed figure shows that e m differs only by approximately 0.39z/ from the 
theoretical limiting value 1/(3-71"). The curve in this figure shows the least-squares quadratic 
fit (1/(3tt) - €m)/u = 0.3911 + 0.9102// + 40.50z/ 2 to the numerical results indicated by the 
diamonds. 

We now discuss the results from a second set of simulations, differing from the first 
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Figure 5: least-squares quadratic fit (1/(3tt) - e m )/v = 0.3911 + 0.9102z/ + 40.50z/ 2 to the 
numerical results indicated by the diamonds. 

only in the initial condition: «2(0) = — 1 and Uk{0) = for k ^ 2. In physical space this 
corresponds to u(x, 0) = — sin2x. For this case, only even wave numbers can be excited. 
Initially, the steepest slope is —2 occurring at x = ±7r, 0, where the inviscid solution blows 
up simultaneously when t — T — 1/2. One would expect e m to be twice as great as that 
in the previous case because the combined contribution to e m at both x = —tt and x = it 
is equivalent to that at x = 0. Furthermore, since the local extrema are 7r/4 away from 
the (stationary) locations of wave breaking, one would expect £* = tt/A. These are actually 
what we have observed. More precisely, the numerics have returned e m = 0.2121 and 
= 0.7856. The spectrum plot is the same as figure [3] and is not shown. 
In passing, it is worth mentioning that for the present example, e m can be made ar- 
bitrarily large by changing the initial condition. Given ue(0) = —1 and Wfc(0) = for 
k ytz £, which corresponds to u(x,0) = — sin £2; in physical space, only the wave numbers 
£, 2£, 3£, ■ ■ ■ can be excited. Initially, the steepest slope is — £ occurring at x = Inn I £ for 
n = 0, ±1, ±2, ■ • ■ and \n\ < £/2, where the inviscid solution blows up simultaneously when 




t = T = l/£. The local extrema are ir/(2£) away from the (stationary) locations of wave 
breaking. One can expect e m = £/{3n) and = ir/(2£), which we have actually observed 
(within small errors as the cases reported above) for several different values of I. Note that 
although e m can be made arbitrarily large by increasing £, Eq. (TH1) does hold as both of its 
sides are proportional to £ (ko = £/L). The scaling E(k,t*) = Ck~ 2 , starting from k = £, 
has been observed to prevail for all cases, with C oc £. 

5 Conclusion 

In summary, we have studied both analytically and numerically one- dimensional viscous 
Burgers flows decaying from smooth initial conditions. The results obtained include upper 
bounds for the energy dissipation rate and number of degrees of freedom and constraints on 
the spectral distribution of energy. Given that the maximally achievable energy dissipation 
rate e m remains finite and positive in the inviscid limit v — > 0, it is found that energy 
spectra steeper than k~ 2 are ruled out in that limit. For this critical scaling, e m satisfies 
e m < V3ko Hiiolloo H M o|| 2 , where ko is the lower wave number end of the energy inertial range 
and u is the initial velocity field. This further implies the upper bound k v < Hwolloo jv 
for the energy dissipation wave number k v . It follows that the number N of Fourier modes 
within the energy inertial range satisfies N < \[?>L HmoIIoo/^? where L is the domain size. 
This result coincides with a rigorous estimate, using no assumption of power-law spectra, 
for the number of degrees of freedom D defined in terms of local Lyapunov exponents. 

As an illustrative example, we have considered both analytically and numerically the 
Burgers equation in the periodic domain [—ir, n] with the initial condition Uq(x) = — sinx. 
In the former approach, we have tightened up the constraint on the spectral distribution 
of energy by pointing out that no power-law energy spectra other than k~ 2 are realizable. 
A detailed examination of the (unique) limiting weak solution has provided an explanation 
why the maximum velocity is better conserved than the energy. In the latter approach, we 



have demonstrated the exact k~ 2 scaling and have numerically determined the viscosity- 
independent dissipation rate and time of maximum energy dissipation. These are consistent 
with analytic results derived from the limiting weak solution. 
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